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Abstract 

Long distance contribution to the — mixing is taken into account con- 
sistently and the corrections to the naive duaUty result represented by the famous 
box diagram are found to be small. Estimates are given in the leading order of the 
chiral perturbation theory. 
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Both experimental and theoretical investigation of the — system is still drawing 
much attention since it can be used for precise tests of the standard model of particle 
interactions and for searching new physics beyond the standard model. Mixing of neutral 
mesons was calculated in |]l| in the leading order of perturbation theory (PT) having 
given also a successful prediction for c-quark mass through GIM mechanism . However 
it has later been pointed out that the local effective AS = 2 Lagrangian cannot account 
for the long distance contribution which is present in the initial Green's function for the 

_ j^o ]2iixing and is connected with the propagation of the light w-quark round the 
loop of the box diagram ^. Within effective theory approach transitions with AS = 2 
are given by the second order of S-matrix expansion in the effective AS* = 1 Lagrangian 
that is given by colorless four-quark operators. Long distance contributions are then 
related to the fact that the corresponding two-point correlator of four-quark operators 
reveals a nonperturbative behavior at large distances due to formation of intermediate 
meson states P, ^ ^, ^, |^. The estimates of the hadronic contribution obtained in this 
way suffer from large uncertainties because at present the relevant matrix elements (or 
form factors) cannot be directly extracted from experimental data with sufficient accuracy 
(see, for example, ref. |^). Moreover a consistent separation of long and short distance 
contributions is not thoroughly defined in the standard approach and hadronic contribu- 
tions are added independently to the result of perturbative calculations. This can lead 
to a double counting because of duality between quark and hadron description of mix- 
ing. Indeed, due to GIM cancelation the main contribution to the PT result comes from 
ti-quark with a momentum smaller than rric running inside the loop of the box diagram. 
On the other hand this perturbative QCD part after inclusion of non-perturbative QCD 
corrections has to reproduce all the hadronic contribution below the scale rric according 
to QCD-hadron duality. Furthermore, since the perturbative result essentially saturates 
the experimentally observed Kl — Ks mass difference one can conclude that the PT 
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contribution is dual to the hadronic one up to small corrections. In the present paper we 
estimate corrections to this naive duality relation. 

The tree level local effective AS = 1 Lagrangian after decoupling the heavy particles 
(ly— boson, t and b quarks) has the form 



G 

Las=i = -^Vsq'V^gS-f^{l - 75)g'g7M(l - 75)c?, g', q = u, c (1) 



where is the Fermi constant, V is the Cabibbo-Kobayashi-Maskawa matrix of quark 
mixing. 

The matrix element M of the K'^ — transition is represented by the expression 

M='-J dx{K\k')\TLAs=i{x)LAs=imK\k)). (2) 
The effective AS* = 2 Lagrangian follows from eqs. 

+ (3) 

Here 9c is the Cabibbo angle that parameterizes the quark mixing in the limit of two 
generations to which we restrict our consideration in the following. The heavy part of the 
transition 

Lh = i J THix)dx, Th = Tec - Tcu - T^c, 
T,,. = (K°|rs7;.(l - l,)qq'l,{l - l,)d{x) S7,(l - l^)q'qiA^ - 75)^(0)1^°), (4) 
contains loops with virtual heavy c-quark in the intermediate state while 

LL = i J TL{x)dx, Tl = Tuu, (5) 

describes the light part of the transition. 

Expression (|^) is finite due to GIM mechanism but Lh and separately require 
regularization because they are ultraviolet divergent. We use dimensional regularization. 



At first sight, an explicit cutoff in the momentum space seems to be more relevant in 
the context of GIM cancelation. However, it violates chiral symmetry which is the basis 
of our analysis of low energy dynamics of strong interactions within chiral perturbation 
theory (ChPT). 

The heavy part of the Lagrangian @ with virtual c-quark can be computed within 
ordinary PT. It corresponds to localization at distances of order 1/mc and the result reads 

LH = -^,{s7,{l-l,)d)' (6) 

that gives the standard contribution to the matrix element (0) 

/ Gf sin cos / l\ml 



where X = f%m\-, 



{s^,i,df = 2(1 + ^) BKmlfl, 



and Bk = 1 in the vacuum saturation approximation. We consider u and d quarks as 
massless ones. 

Corrections to this result are under control. The leading PT corrections due to strong 
interactions has been obtained long ago in ref . p[ . Recently the next-to-leading perturba- 



tive calculation has been completed f^, |T0[- The hadronic matrix element of the four-quark 
operator (D has been intensively studied during several last years with different techniques 
|TT| and the final result is consistent with the vacuum saturation hypothesis within a 20% 



accuracy. Regular corrections in l/iric i.e. the contribution of operators with dimension 
eight in mass units have also been taken into account and happened to be small |T^. 

Turning to the light part of the amplitude (^), it is not represented by the PT expansion 
because long distance effects have to be taken into account. If one assumes the validity of 
formal PT expansion and neglects the nonperturbative effects of {K^) state formation 
the result for this part vanishes in the chiral limit within dimensional regularization. The 
aim of the present paper is to calculate the nonperturbative deviation from this result. 



To estimate the light part of the amphtude we use the factorization procedure, i.e. 
we suppose the K'^ and K'^ states to be interpolated by the two-quark operators with the 
wave function Wfj,{k,x) [|l3l 



The factorization approximation is used in the standard analysis of the box diagram 
contribution and we are interested rather in the corrections connected with the non- 
perturbative behavior of the w-quark propagator at large distance than in the violation 
of this hypothesis. Adopting this approximation one has a wave function of the kaon 
factored out even for separated space-time points. The light part of the amplitude can 
now be written in the following form 



Two terms in the right hand side of this equation are the direct and cross parts respectively. 
In the direct contribution a pair of s and d quarks which interpolate {K^) state is 
localized in a space-time point and we have 
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The wave function in eq. (P) is normalized according to the equation 






Ik^k.i I dxe'''^{0\Tuj^{l-^,)u{x)u^,{l-^,)umO). (9) 



Using the anomaly relation d^{u'~^^'~^^u) = 2w where 



w{x) 



1 
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eq. can be rewritten in the following form 
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In the small k limit an expansion of the two-point correlator of w operators is known 
within the chiral perturbation theory (ChPT) 

z I dxe''^Q\Twix)wiO)\0)\,.^^2^ = \ml{fl + ^o) + O(m^) 

where Hq is the vacuum susceptibility 



We have used a mass relation 



which is valid for massless u and d quarks. Finally the direct contribution becomes 

2/1^2 



= 3(^) W.' + ^o). (10) 

In the cross contribution s and d quarks which interpolate state come from 

different space-time points and we have 

where 11* and 11° are singlet and octet contributions respectively 

= ^1 cix(0|TM7,(l - 75)m(x)u7^(1 - 75)^(0) |0)(2^/,„(7,;3 - g^,ugo.(i)Wl{k,x)W,{k,x), 

IT = i J rfx(0|TM7at''(l - 75)M(x)M7/3t''(l - -f5)u{0)){2gf,a9uf3 - gf,u9ai3)W*{k, x)W„{k, x), 

and t"- are SU{3)c generators normalized with the condition tr{t°'t'') = 6"'^/2. If the PT 
expression is used for u-quark propagators this contribution computed in the dimensional 
regularization is non-zero even in the chiral limit because the kaon wave function intro- 
duces the nonperturbative scale ~ ^qcd- Note that this contribution is suppressed by 
the factor {^Qco/fTicf' in comparison with the leading one and is not taken into account 
in the standard analysis. Clearly, an estimate of the contribution of massive states to the 
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cross part of the amplitude depends on the model for the ii'-meson wave function. At 
the same time the contribution of the vr^ and t] mesons to the color singlet part 11'* in the 
chiral limit is proportional to \Wfj,{k, 0)p and is model independent 

Contribution of states with nonzero masses are suppressed by the ratio of the size of the 
wave function to the mass of the state and are neglected. One can reasonably suppose that 
the magnitude of long distance effects in this channel is semiquantitatively represented by 
this simple estimate. The advantage is also the additional suppression of the nondirect 
contribution in A^^ counting (it is suppressed by the factor 1/Nj}) that allows one not to be 
very demanding to the accuracy of this estimate. In fact, we assume that the contribution 



of higher mass states can amount about fifty percent of the estimate (|T^) still not leading 
to the large uncertainty of the total result. 

Note that if the explicit cutoff in momentum space /i is used in calculation of the light 
part of the amplitude for cutting high energy contributions, the heavy part (Q) of the 
full matrix element has to be properly modified. For fi = rric the heavy part vanishes and 
the PT result for the light part coincides with eq. (|^). On the other hand the light part of 
the amplitude can now be computed by saturating eq. @ by contributions of the hadrons 
with masses smaller than rric what is the standard way to estimate the long distance 
contribution. However in this case the light part (or long distance contribution) is not 
an addition to (|^) but rather represents the whole amplitude. So the nonperturbative 
corrections in this case should be associated with the difference between the hadronic and 
PT expressions for the light part of the amplitude not with the hadronic contribution 
itself. 

We would like to stress this point. The physical result being invariant, its repre- 
sentation through (and splitting into) long and short distance contributions depends on 
regularization procedures used for computation Lh and L^. In these circumstances the 



advantage of dimensional regularization for perturbative computation of the light part of 
the whole amplitude becomes clear - the light part vanishes. The dominant contribution 
to the amplitude stems from the heavy part which can be reliably computed within PT 
while the light part contributes only because of some nonperturbative effects. These non- 
perturbative effects under an assumption of factorization are organized into a correlator 
of colorless currents of light quarks taken at the scale of order m|-. ChPT is now used to 
compute the final answer that is therefore exact. Numerically it is rather small. An ad- 
ditional advantage that makes the result more reliable is that the correlator is suppressed 
in formal counting in Nc- If, instead, the explicit cutoff in momentum space is used for 
regularization, hadronic contributions represent the total result. Because the energy scale 
rric is rather large for light quarks the spectral density cannot be accurately given that 
leads to large uncertainties 0. So, it seems impossible to estimate the corrections to 
eq. (0) in this way. Moreover the multihadron contributions generate a power dependence 
on the cutoff (see, for example, ref. 0) and for ^ ~ there is no reliable way to find 
their numerical value. On the other hand if the cutoff is small and chosen to be of order 
of the chiral symmetry breaking scale fi ~ itlk, the contribution of hadronic states from 
momentum region rrix < k < rric is out of control because it can be described neither 
within pQCD no ChPT. This is the reason for considerable uncertainties of previous es- 
timates of the hadronic contribution to the mixing. At the same time in the dimensional 
regularization we deal with the well defined chiral expansion in m|^/A^ ~ 0.25 where 
A^ ~ 1 GeV is the scale of ChPT. 

The obtained results are applied to calculation of the parameters of — system. 
The correction to the Kl — Ks mass difference reads 

8^V^ Ho I W ^ f-, 0.03\ 

where Am^rf is the short distance contribution that follows from eq. and the correction 
comes from eqs. ( ]T0| , |I2|). For numerical estimate we use = 1.4 GeV, /^r = 0.132 GeV. 



Thus we found that the difference between the full result for the matrix element of 
the mixing and the contribution of the heavy part computed within the dimensional 
regularization is about 3% of the total result. This difference vanishes in large Nc limit. 

To conclude we would like to stress again the main result of the paper. We showed that 
the contribution of light hadrons to K'^ — transition amplitude (known as dispersive or 
long distance contribution as well) should not be added independently but is essentially 
dual to PT expression of the box diagram with virtual xi-quark. In our analysis we rely 
on factorization approximation even for nonlocal operators and intensively use the chiral 
limit of QCD. Within these assumptions duality is rather accurate quantitatively and 
the PT expression represents the complete answer well while the genuine nonperturbative 
contribution of light quarks to mixing is negligible (about 3%). It has been proven that 
the same result can be also obtained by direct summation of hadronic contributions till 
energies of order rric- The latter method is however impractical at present because of poor 
experimental data. 
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